No quasi-long-range order in the two-dimensional liquid crystal 
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Systems with global symmetry group 0(2) experience topological transition in the 2-dimensional 
space. But there is controversy about such a transition for systems with global symmetry group 
0(3). In this paper, we study the Lebwohl-Lasher model for the two-dimensional liquid crystal, using 
three different methods independent of the proper values of possible critical exponents. Namely, we 
analyze the at-equilibrium order parameter distribution function with: I) the hyperscaling relation; 
2) the first scaling collapse for the probability distribution function; and 3) the Binder's cumulant. 
We give strong evidences for definite lack of a line of critical points at low temperatures in the 
Lebwohl-Lasher model, contrary to conclusions of a number of previous numerical studies. 
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Introduction. Mermin and Wagner [l[ established that 
no ferromagnetic phase nor any long range order can ap- 
pear for systems of continuous symmetry at finite tem- 
perature in space dimension d < 2. However, such sys- 
tems might have another type of transition governed by 
binding-unbindingtopological defects at definite positive 
temperature T BKT 0, y, 0| ■ This kind of topological phase 
transition is called Berezinskii, Kosterlitz and Thouless 
(BKT) transition. 

The two-dimensional (2d) XY-model, with global sym- 
metry group 0(2), exhibits such topological transition (ij. 
Quasi-long-range order (QLRO) appears at low temper- 
atures T, and the order parameter vanishes as a power 
law at the thermodynamic limit. Due to the QLRO be- 
havior, the system susceptibility %, that measures the 
fluctuations of the order parameter, diverges for all tem- 
peratures T < T BKT , and the system is characterized by 
a line of critical points below the critical temperature 
T BKT . Close to T = 0, correlations are dominated by 
spin-wave solution: in the units system where k B = 1 
and the coupling factor between magnetic moments is 
J = 1, the correlation function exponent, 77, depends on 
the temperature as: 77 = T/2ir. Another characteristic 
behavior of this transition is that at temperatures just 
above the BKT transition, t = (T - T BKT ) /T BKT > 0, the 
correlation length, £, diverges as the essential singularity: 
£ ~ exp(6< -1 / 2 ), that is much strongly than the ordinary 
second order transition power law, £ ~ t~ v '. 

On the other hand, Polyakov [B[ , using renormalization 
group theory, proved that the 2d Heisenberg model, with 
global symmetry group 0(3), does not present any sort 
of phase transition. An important difference between the 
two systems above is that the global symmetry group for 
the XY-model is abelian, while it is not for the Heisen- 
berg model. However, things are not as clear: numerical 
evidences were recently given for transition S, Sfi @, E3] 
in this system, and a possible QLRO phase [y, at very 
low temperatures. In the same spirit, it has also been 



reported that the 2d fully frustrated anti-ferromagnetic 
Heisenberg model presents a crossover produced by the 
binding-unbinding of topological defects in a very nar- 
row temperature interval. In this case no QLRO behav- 
ior at or below the transition I^, H| has been observed. 
At present, this controversy on systems with continuous 
non-abelian symmetry is not solved. 

Kunz and Zumbach (KZ) [3] performed intensive 
Monte Carlo study of the 2d RP 2 model, which has 
the global symmetry group 0(3) and the local symme- 
try group Z2- The model describes the isotropic- nematic 
transition of a liquid crystal. KZ concluded with BKT- 
like transition from analysis of energy, specific heat and 
topological quantities. But the correlation length behav- 
ior for t > was not proven to be either of the power-law 
or of the essential singularity type. Ten years later, the 
problem of phase transitions for liquid crystals in d = 2 
was complemented [1, [TjJ [lj| using the powerful tech- 
niques of conformal transformations [HI (CT) and finite 
size scaling (FSS) [EGl]. The Lebwohl-Lasher (LL) 
was preferred to the RP 2 model though sharing the same 
symmetries. These studies concluded with a BKT-likc 
transition and a QLRO phase below the BKT tempera- 
ture estimated by KZ 14[ to be T BKT = 0.513. At low 
temperatures, a spin wave dependence 77 oc T was ob- 
tained. To support this conclusion Dutta and Roy [l9l | 
showed that the transition is driven by topological stable 
points defects known as i-disclination points. Using the 
system susceptibility, x, FSS was able to estimate the 
value of the correlation function exponent 7/ within the 
temperature range T < T BKT . On a line of critical points, 
X should scale with the exponent 7/f, which is related 
to 77 through the hyperscaling law: 



7/1/ = 2 — 77. 



(1) 



Us ing, ([T]), estimation of the values of 77 was performed 
[(J [20] ■ The values appeared to behave similarly to the 
ones obtained through CT and scaling of the order pa- 
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rameter, but there is a discrepancy of about 5% between 
both results. The origin of such difference was tentatively 
explained arguing that the system sizes were far from the 
thermodynamic limit and the number of independent re- 
alizations were too small to reach good statistics. 

The purpose of this article is to revisit the problem 
of the possible appearance of BKT-like transition for the 
2d LL model. In this model, liquid-crystal molecules are 
represented by unitary vectors <7j situated on the sites, 
labelled i, of a hypercubic lattice A of length L. The 
Hamiltonian is given by: 



We shall use now alternative method to check the hy- 
perscaling relation. Let us introduce a as the standard 



0H = ^^P 2 (5W* 



+s) 



(2) 



where (3 — 1/T, Pi is the second Legendre polynomial 
and the interaction is between nearest neighbors. The 
appearance of the P2 function in @ comes from the Z-i 
local symmetry. In the nematic phase, the preferential 
direction is characterized by the unit vector n, called the 
director, and one can measure the local orientation with 
respect to the director by: <jj • n = cos^. Then, the 
local order parameter is defined by m(i) — ^(cos^j)}. 
Whenever the system is completely ordered: m(i) = 1. 

Therefore, we start from the hypothesis that the 2d 
LL model experiments BKT transition at T BKT = 0.513, 
similarly to the transition observed in the 2d XY-model 
at T BKT = 0.893. Below such critical temperature a line 
of critical points should be observed in both models. 
To validate this point, we performed Monte Carlo 
simulations using the Wolff algorithm 21] in d = 2, with 
periodic boundary conditions at temperatures well below 
T BKT . A total of 6 x 10 6 independent realizations were 
performed for each system size and each temperature 
for both models. Then we found estimates of the order 
parameter probability distribution function (PDF), and 
estimation of the validity for the hyperscaling relation 
(jTJ) . Finally we will analyze the Binder's cumulant 
behavior, comparing also with the Heisenberg model. 
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deviation of the order parameter. One has 
then a scales with the system size as: a ~ 
the 2d systems. Therefore, the ratio (m)/a should be 
a constant whenever the hyperscaling relation (jTJ) is sat- 
isfied. The great advantage for using this ratio is that 
previous estimation of the exponents is not necessary to 
check (JT). In Fig. [TJ {m)/a is plotted versus L 1 for 
the XY (above) and the LL (below) models in the low 
temperature domain. 

For the XY-model, the ratio is seen to saturate at the 
thermodynamic limit to a value (m)/o~ ~ 31.1. 

For the LL-model, power law is the best fit consistent 
with our data for (to) /a versus L. The ratio does not 
saturate to a finite value and we conclude that the hy- 
perscaling relation (fT]) does not hold in this case. 

Similar behavior was observed for the XY-model at 
TUt = 0.893 [IH and for the LL-model at T BKT = 0.513 
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Hyperscaling relation check. - For the XY-model at 
T = 0.6, we observe that both the order parameter and 
the susceptibility have power law behavior, (to) ~ L~PI V 
and x ~ L 1 ^ respectively. For the XY system the expo- 
nents obtained were (3/v = rj/2 « 0.058 and 7/^ m 1.877. 
With use of the CT method, Berche et al [22| obtained 
the value (3/v = 0.0595 in excellent agreement with our 
results. Hyperscaling relation (TTJ is satisfied with error 
smaller than 0.4%. 

For the LL model at T = 0.4, we obtained again excel- 
lent power laws for (to) and x with respective exponents 
P/v = 77/2 w 0.0945 and j/u « 1.868. But now, the 
agreement for Eq.(fT]) is poor and about 3% (one order 
of magnitude larger than for the XY case). The actual 
increase of the number of independent realizations does 
not really improve the results obtained previously @, |20| . 



FIG. 1: {color online) {m)/a is plotted vs L 1 for the XY- 
model at T = 0.6 (top) and for the LL-model at T = 0.4 
(bottom). A linear fit is obtained for the XY-model ((m) ja = 
31.1 — 16.4/L). A power law fit shows that no saturation 
is observed for the LL-model. Both fits are shown as bold 
lines. The bold circles are the data from The number 
of independent realizations used to obtain the bold squares 
is almost two orders of magnitude larger than in [(|. The 
hyperscaling relation (jTJ is not satisfied in the thermodynamic 
limit by the LL-model at T = 0.4. 

First-scaling relation check. - The first-scaling law 

Tfl 

(m)P(m) = 4>t(zi), with z\ = , (3) 

< to > 

should be satisfied anywhere on the line of critical points 
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below the BKT transition. In Q, P(m) denotes the or- 
der parameter PDF. The scaling function $t depends 
only on the actual temperature. Here too, one great ad- 
vantage of the first-scaling law is that Eq.© does not 
require knowledge of any critical exponent. In Fig. [2] the 
order parameter PDF is plotted in the first-scaling form 
for both models. 

For the XY-model, the three curves exhibit almost per- 
fect collapse. Relation ([3]) is clearly satisfied at T = 0.6. 
Similar behavior was observed previously for the XY- 
model at T BKT [23j . The definite shape of the scaled dis- 
tribution is Weibull-like [2(| similarly to the T BKT case 

M- 

For the LL model, collapse is not realized in Fig. [2] 
As the system size is increased, the scaled distributions 
tend to separate for T = 0.4. This is evidence that the 
LL model is not at a critical point for this temperature. 



is known to be a universal quantity independent on L at 
the critical point 25 1. 

For the XY-model, t/ 4 is universal for T < T BKT [2f|. It 
is checked on the Fig. [3]where U4 is plotted for this model 
(above). For the XY-model a crossing point is observed 
near the reported BKT temperature. For temperatures 
below the crossing point, the U4 grows with the system 
size. All the curves are expected to collapse in this inter- 
val when L — > 00. It is faster when the temperature is 
small [2(|. On the other hand, the U4 above the crossing 
point, decreases with increasing L. This type of behav- 
ior is observed in others 0(2) models with Z2 symmetry 
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The behavior of U4 is completely different for the LL 
model. The Binder cumulant decreases with L in all the 
domain of temperature explored (T > 0.1). No crossing 
is observed anywhere in Fig. [3] (bottom). 
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FIG. 2: (color online) Order parameter PDF for the XY- 
model at T = 0.6 (top) and the LL model at T = 0.4 (bottom) 
in the first scaling form. A perfect collapse is observed for the 
XY-model. This is not the case for the LL model. The L = 
768 data are from [jj, with 9 x 10 4 independent realizations. 
It is clear from this figure that the number of independent 
realizations used in [6J was large enough to realize the first- 
scaling law. No self-similarity is observed for the LL model 
at T = 0.4. 

Binder's cumulant check. - For a continuous phase 
transition the Binder's cumulant, 
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FIG. 3: The Binder cumulant vs the temperature for the XY- 
model T — 0.6 (top) and the LL model (bottom). The number 
of independent realizations is 10 5 for each T and L. No cross- 
ing is observed for the LL model. Therefore, no evidence of 
any phase transition is observed for the LL model. Values of 
the Binder cumulant for T = 0.4 are shown on the inset as 
a function of the system size. We used 6 x 10 6 independent 
realizations to obtain each point, so that the error bars are 
much smaller than the symbol size. 

To complement the discussion, we study the Binder 
cumulant behavior for the Heisenberg model at T > 0.1. 
It is seen to behave very similarly to the 2d LL model, 
as no crossing is observed (see Fig. @|. Then we can con- 
clude that in the low temperature range, the LL model 
must have very large (but not infinite) correlation length, 
that suddenly begin to decrease in the neighborhood of 



4 



T = 0.513. For this reason an apparent QLRO phase may 
be observed below this temperature. 
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FIG. 4: The 2d Binder cumulant for the Heisenberg model 
exhibits the same type of behavior as for the LL model (see 
Fig. |3J. 

Discussion. - We presented in this paper three strong 
evidences supporting the idea that the 2d liquid crystals 
do not have a quasi-long-range order phase, namely: 

(a) the hyperscaling |T]) is not satisfied; 

(b) the first-scaling collapse J3]) does not hold; 

(c) the Binder cumulant Q does not exhibit any crossing 
point. 

Then this system can not experience a transition of the 
BKT type. 

From FSS analysis, Mondal and Roy [29( concluded 
that the LL model should present a continuous transi- 
tion at T = 0.548. The lack of crossing event for the 
Binder cumulant behavior (as observed in Fig. [3]) def- 
initely suggests that this is not the case. In reference 
0, Ho] the stiffness and the susceptibility are studied as 
functions of temperature T and system size L for the 
LL- and the XY-models. For the XY-model the stiffness 
saturates to finite value below T BKT . However for the 
LL model the stiffness tends to decrease logarithmically 
with the system size, similar to the behavior of the fully 
frustrated anti-ferromagnetic Heisenberg model (FFAH) 
[l3j |. On the other hand the susceptibility for the LL 
model changes its functional form in a small region of 
temperature around T = 0.513. This is also observed in 
the FFAH [13{ . Then for this reason, and knowing the 
fact that topological defects are stable [l9[ , we speculate 
that the LL model may have a crossover similar to FFAH. 

The set of critical- exponents free methods used in this 
article can be used to explore any thermodynamic sys- 
tems and to identify possible critical points. The hy- 
perscaling relation and the first-scaling law are of great 
utility to identify whether a system is or is not at a criti- 



cal point. In particular, such procedure could be helpful 
for the Heisenberg model in the T < 0.1 domain, to dis- 
cuss about a p ossible transition at very low temperature 

AIFS and RP like to thank Bertrand Berche for stim- 
ulating discussions about the subject of this article. 
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